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QUASI-STATIC  RANGE  PROPAGATION  EQUATIONS 
FOR  THE  APPROXIMATE  FIELDS  WITHIN 
A CONDUCTING  SlAB 


PRODUCTION 

During  the  last  several  years,  considerable  Interest  has  developed  in 
detennining  the  quasi-static  field  components  of  antennas  located  above,  or 
buried  within,  the  earth's  surface.  Quasi-static  range  is  defined  as  the 
range  of  transmission  where  the  measurement  distance  is  much  less  than  a 
free-space  wavelength.  Quasi-static  range  results  are  useful  for  submarine 
radio  conmtunication  and  detection;  they  also  have  applications  in  locating 
buried  miners  and  help  geophysicists  determine  the  electrical  properties  of 
the  earth. 

For  the  semi-infinite  conducting  medium  case  (i.e.,  air,  single  layered 
earth),  some  work  has  been  done  to  determine  the  quasi-static  fields  pr^uced 
by  various  subsurface  sources  when  the  measurement  distance  is  comparable  to 
the  earth  skin  depth.'"*  (Most  of  these  results  are  summarized  by  Kraichman.^) 
However,  the  field  strength  expressions  are  very  coaplex  because  they  involve 
products  of  modified  Ressel  functions  of  different  argument.  Recently,  by 
using  finitely  conducting  earth-image  theory  techniques,  we  derived  approxi- 
mate expressions  for  the  general  quasi-static  range  electromagnetic  fields 
(in  air  and  in  earth)  produced  by  various  subsurface  antennas.®  Some  numer- 
ical calculations  have  also  been  provided.® 

These  investigations  revealed  that  when  the  source  is  not  buried  too  far 
from  the  surface  of  the  conducting  half-space,  the  resultant  electromagnetic 
field  is  significantly  different  from  that  obtained  in  an  infinite  conducting 
medium.  Therefore,  if  the  conductor  is  two-layered,  the  lower  Interface  must 
likewise  affect  the  field  (especially  if  the  upper  layer  in  which  the  source 
is  located  is  not  too  deep).  For  example,  if  the  transmitting  antenna  is 
situated  in  a shallow  sea,  the  theory  for  a uniform  conducting  half-space 
would  not  accurately  describe  the  field  at  the  sea  bottom.  This  observation 
is  also  true  at  the  surface  of  the  sea  if  the  sea  depth  (Ij)  is  less  than 
approximately  one  skin  depth. 

Weaver  has  evaluated  the  exact  Sonmierfeld  integral  expressions  and 
obtained  numerical  results  for  the  quasi-static  fields  produced  by  horizontal 
and  vertical  electric  dipole  (HED  and  VED)  antennas  located  in  the  upper 
layer  of  a two- layer  conducting  half-space.®  Numerical  results  for  the  quasi- 
static fields  produced  in  the  sea  by  a vertical  magnetic  dipole  (VMD)  for 
various  values  of  sea  bottom  conductivity  have  been  obtained  by  Coggon  and 
Morrison. 

This  report  employs  finitely  conducting  earth-image  theory  techniques  to 
derive  approximate  expressions  for  the  general  ac  quasi-static  range  electro- 
magnetic Helds  produced  by  various  subsurface  dipole  antennas,  which  include 
the  HED,  VED,  VMD,  and  horizontal  magnetic  dipole  (HMD)  types,  each  of  which 
is  located  in  the  upper  layer  of  a two-layer  conducting  earth.  For  mathema- 
tical convenience,  consider  the  conductivity  of  the  bottom  layer  as  being 
eqttal  to  zero.  Thus,  the  problem  is  reduced  to  calculating  the  fields  in  a 
conduct ing slab.  Treat ing  the  bottom  layer  conductivity  as  eqxial  to  zero  is  not 
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that  restrictive; in  many  practical  cases,  the  conductivity  of  the  upper  layer 
(oi)  is  much  greater  than  the  conductivity  of  the  bottom  layer  (02),  particu- 
larly in  the  sea-to-sea  bed  case.  Therefore,  this  assumption  limits  the 
results  to  measurement  distances  (p)  of  approximately  62/-‘’«  wh®re  ^2 
skin  depth  in  the  bottom  layer.  For  example,  at  a frequency  of  1 Hz,  if  oj 
* 4 S/m  (sea)  and  02  ■ 10'^  S/m  (sea  bed),  then  6j  "v  250  m and  62  5 km. 

Therefore,  for  this  example,  the  results  should  be  valid  to  a measurement 
distance  of  approximately  1 km.  Furthermore,  if  the  conductivity  of  the 
sea  bed  is  10“^  S/m  (62  50  km),  the  results  should  be  valid  to  a measurenent 

distance  of  approximately  10  km. 

For  the  purpose  of  this  report,  all  four  sources  are  located  at  depth  h 
(h  i 0)  with  respect  to  a cylindrical  coordinate  system  (p,^,z)  and  are 
assumed  to  carry  a constant  current,  I.  The  VED  and  HED  antennas  (of  infin- 
itesimal length  t)  are  oriented  in  the  z and  x directions,  respectively.  The 
axes  of  the  Vk®  and  HMD  antennas  (of  infinitesimal  area  A)  are  oriented  in  the 
I and  y directions,  respectively.  Free  space  occupies  the  regions  z < 0 and 
z > ti,  whereas  the  conducting  slab  occupies  the  region  0 < z < ti(see  figure 
1).*  Displacement  currents  are  neglected  in  both  the  slab  and  the  air.  The 
magnetic  permeability  of  the  conducting  slab  is  assumed  to  equal  pg,  the 
permeability  of  free  space.  Meter-ki logram-second  (mks)  units  are  employed 
and  a suppressed  time  factor  of  exp  (ioit)  is  assumed. 


MODIFIED  FINITELY-CONDUCTING, 

EARTH- I MAGE  THEORY  TECHNIQUES 

In  the  semi-infinite  conducting  medium,  where  both  the  source  and  receiving 
antennas  are  located  above  the  earth's  surface  (z,h  < 0)  (see  figure  1),  the 
quasi-static  range  (yg  0)  integrals  are  of  the  type  that  can  not  be  evalu- 
ated analytically  throughout  the  quasi-static  range,**  that  is 


/ (^) 


e J^(Xp)dx. 


(1) 


where 


Y 

u 


argument 


= iu)(ppep)^'*^^  'V'  0 (air), 

1/2 

“v  (iuppo)  (earth) , 

. 2 

■ (X  ♦ Y ) , and 

j (Xp)  * Bessel  function  of  the  first  kind,  order  zero,  and 
Xp? 


•Figures  1 through  6 are  presented  at  the  end  of  the  text. 
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Physically,  the  essence  of  the  finitely  conducting  earth-image  theory 
technique  is  to  replace  the  finitely  conducting  earth  with  a perfectly  con- 
ducting earth  located  at  the  (complex)  depth  A/l,  where  d » 2/y  - 6(l-i). 
Analytically,  this  corresponds  to  replacing  the  algebraic  reflection  coeffi 
cient  (u-X)/(U'»X)  in  the  exact  integral  equations  by  exp  (-Xd),  where  X is 
the  variable  of  integration.  Once  this  is  accomplished,  (1)  can  be  readily 
evaluated.  For  antennas  located  at,  or  above,  the  earth's  surface,  the 
general  image  theory  approximation  is  valid  throughout  the  quasi-static 
r.utge.^ 

When  the  source  and  receiving  antennas  are  located  below  the  earth's 
surface  (z,  h > 0)  (see  figure  1),  the  quasi-static  range  integrals  to  be 
evaluated  are  of  the  form*'®'® 


' / (Hr)  » 


■u(r  ♦ h) 


J {:ip)dx 


(2) 


Integrals  of  this  type  can,  and  have,  been  evaluated  analytically 
throughout  the  quasi-static  range.**®  However,  the  resulting  expressions  are 
very  complex  because  they  involve  products  of  modified  Bessel  functions  of 
different  argument.  Therefore,  wc  let® 


-u(z  h)  -Ya(z  h)  -Yb(z  ♦ h) 
e "v  e e 


(3) 


where 

a ■ 0 and  b » 1 for  Rj/6  <<  1 |Rj  “ ♦ (z  + h)*”  j , 

a - 0.4  and  b = 0.96  for  R^/fi  less  than  approximately  1, 

a « 0.96  and  b •=  0.4  for  R^/S  between  approximately  1 and  10,  and 

a ■ 1 and  b » 0 for  p > 3(Z''h). 


Substituting  (3)  into  (2)  and  substituting  exp(-Xd)  for  (u-X)/(u+X) 
results  in 

^ . h)  J ,-x|d  . b(z  . ^ 

o 

yhich  can  be  readily  evaluated  throughout  the  quasi-static  range. 


DERIVATION  OF  THE  CONDUCTING  SlAB  FIELD- 
COMPONENT  EXPRESSIONS  FOR  THE 
GENERAL  QUASI -STATIC  RANGE 

Because  we  have  already  derived  the  general  quasi-static  range  sub- 
surf acc-to-subsurface  field-component  expressions  for  the  semi-infinite 
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medium  case,  we  can  employ  these  results  and  the  method  of  images  to  derive 
approximate  expressions  for  the  ac  quasi-static  fields  within  a conducting 
slab.®  In  this  situation,  the  multiple  image  pattern  is  an  infinite  array 
that  supplements  the  original  dipole  exactly  as  the  images  of  a physical 
object  located  between  two  plane  mirrors  appear  to  be  an  infinite  array  of 
that  object. 


HORIZONTAL  ELPCTRIC  PirOLL  (HEP) 


The  HEP  equations  can  easily  be  derived  from  equations  (64)  through  (69) 
of  an  earlier  report.®  They  are 
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T.  • ^ “ \ d ♦ bX  bX  -YaX 

„ . . I .1  -5—1  - -5-1-  . ' 

» "■»  ■'  *2N  ' *1N 


4-  ( 1 ♦ >R^„  \ 

V ON  / 


(d  ♦ bX^) 


It  cos 
4it 


» A 

-*  .1  % i-  ( 


* • ’“on  ' 


X -YR,„  d 4 bX,  bX  ( 

V IN  ) 2 1 

•^IN  ml  j I 


V'"Kr 


2 1 1 
3 ■ 3 

^IN  ^2N 


where 


Xi  - 2nt^  ♦ 2 - h. 


^2  * 2ntj  ♦ r ♦ h. 
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2 2 2 
Rl,  - P . X^. 


2 2 2 2 
Ajn-p  -bX^. 


•*•  C*!  ♦ bX^)^,  and 


Eq  » 1,  * 2,  n = 1,  2,  3, 


HORIZONTAL  MAOIETIC  DIPOLE  (HMD) 


The  HMD  equations  can  be  derived  from  equations  (76)  through  (81)  of  an 
earlier  report.®  They  are 


i(jou„IA  cos4i  « 


E - 
P 


L "n  n3  ( 


n=0  / R 


'"""on  ® 


^2  / X ‘""tn  e’^^2  d > bX  bX 


iuu  lA  sin(Ji  “ 1 X , -YR 

E _ 0 y . /w  YR  ) e ON 

I UN 


(d  + bX^)  -yax^  g‘^®^2  d + bx^  bX^ 
A2n  P^  - ^2N  ^1N_ 


-YR 

iuy„lA  cosiji  * ..ON 

E -v  —S~ I p e 

z 4it  n „3 
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(*  * ’'"on) 


('•''"in)  • 


VERTICAL  ELECTRIC  DIPOLE  (VED) 

The  VED  equations  can  be  derived  from  equations  (73)  through  (75)  of  an 
earlier  report.®  They  are 


E ''' 
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UXDRESSIONS  FOR  TIIF.  QllASI-NFAR  RANGF. 


(22) 


Quasi-ncar  range  Is  defined  us  the  asymptotic  part  of  the  quusi-statlc 
range,  i.e.,  where  the  measurement  distance  (p)  is  much  greater  t!>an  a skin 
depth  in  the  conducting  medium  and  much  greater  than  the  depth  of  burial  of 
the  transmitting  and  receiving  antennas.  Generally,  p must  be  gi'eater  than 
.Xf  and  >3(:  + h) . For  the  conducting  slab  case,  p should  also  be  greater 
than  However,  as  we  shall  see  later,  the  requiix'ment  that  p ^ 23  and 

p 2tj  may  be  sufficient  for  most  cases. 

The  quasi-near  range  approximation  is  setting  the  function  u ■ 
in  the  exact  integi*al  expressions  equal  to  y,  which  is  the  propagation  con- 
stant in  the  conducting  medium. 


LONG  HORIZONTAL  LINE  SOURCIE/ANTFNNA 


Fot'  a long  horiionta 
the  electric  field  within 
I > h) 


where 


u 


line  source  antenna, 
a conducting  slab  may 


oriented  in  the  x-directlon, 
be  written  exactly  as  (for 

Xy  dx,  (2.3) 
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u - u 

R ■ — T — ^ • and 
a u ♦ u 

• o 


* ^0  ^ ^ quasi-static  range. 


The  magnetic  fields  within  the  conducting  slab  may  be  determined  from 

1 1 ^^x 

H « - ^ and  H - t • (26) 

Y luvi  dz  Z iu>p  3v 

o o 

Application  of  the  quasi-static  (yq  '''  0)  and  quasi -near  (u  y)  approxi- 
mations to  (24)  and  (25)  results  in 

. 2 cosh  Y^  — sinh  Yh  I I cosh  Y (*•,  - 2)  ♦ — sinh  Y(f,  - t) 

iHLil  . _L 1 J_L — i 1 1 J-(27) 

u ^ L r,  . 2X  _ 1\ 


I^Y  sinh  Ylj  ^ ^^^ij) 


If  the  slab  is  not  too  thin,  the  denominator  of  (27)  may  be  approximated  as 


Y sinh  yJ^ 


J [l  * ^ coth  r.j]  ' ’ 'S 


Substituting  (27)  and  (28)  into  (23)  and  (26),  and  noting  that* ^ 

00 

/ cos  xy  dx  ■ 0 


results  in 


” 1 

I X cos  Xy  dx  « j 

o y 

K - “S-l  |q^|  (A(z,h)|  , 
itoy 

lly  ^1  |q|  |B(l.h)|  , 
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1 f h ■ 0 and  z ■ Ij * 


A(z,h)  '''  2B(z,h) 


cosh  Ytj  ’ 


(45) 


which  is  identical  to  von  Aulock's  result.*®  (See  chapter  4,  p.  14  of 
Kraichman^  and  Bannister.*®)  Furthermore,  if  h ■ 0 and  tj/iS  > 2,  then 


I 

A(z,h)  2e  |2  cosh  y(*'2  " 

B(2,h)  2e  ^ jcosh  y('*2  - z)  - 2 sinh  y(<*j  ~ ^)  | * 


and 


(44) 


(45) 


which  is  identical  to  von  Aulock's  result.*®  (Also,  see  chapter  4,  p.  1.5  in 
Kmichman.^) 


HORIZONTAl,  tlf.CTRlC  DIFOUi  (HF.P)  ANTUNNA 


By  following  the  same  procedure  outlined  in  the  derivation  of  the  equa- 
tions for  the  long  horizontal  line  source  antenna,  the  quasi-near  range  HED 
antenna  field-component  expressions,  which  are  valid  within  the  conducting 
slab,  may  readily  he  determined.  They  are 


E '\- 
P 


H -v 
P 


and 


H cos  ^ 

2ffop' 

A(z,h)  , 

(46) 

It  sin  Ip 
nop'^ 

A(z,h)  , 

(47) 

0 , 

(48) 

11  sin  ♦ 

3 

ITYP 

Q B(z.h)  , 

(49) 

It  cos 
2syP^ 

i Q B(2,h)  , 

(50) 

3It  sin  ^ 

■ 7 

• A(z,h)  . 

(51) 

2itY  p 


When  h - 0 and  z ■ Ij.  th®  IH-D  antenna  expressions  reduce  to  the  results 
previously  derived  by  Bannister.*^ 


12 


.ii 


TR  5807 


HORIZONTAL  MAGNETIC  DIPOLE  (HMD)  ANTENNA 

By  employing  Maxwell's  equations  and  the  reciprocity  theorem,  we  see 
that  the  quasi-near  range  HMD  antenna  field  component  expressions,  which  are 
valid  within  the  conducting  slab,  may  readily  he  determined  from  the  HED 
antenna  expressions.  They  are 


iuvi  lA  cos  ^ 

E 'V'  2 Q B(h,z)  , 

2.yP^ 


(52) 


iutp  lA  sin  ()i 


E -v 


Q B(h,z)  , 


ITYP 


(53) 


E 0 , 
z 


(54) 


H jA  sin  ^ Q 9B(h,z) 
p 3 Y 3z  ’ 

up 


(55) 


and 


H . lA  -cos  -i  g i_B(h^  zJ 


(56) 


H Q B(h,z)  , 


2nYp 


(57) 


where  B(h,z)  ■ B(z,h),  with  h and  z interchanged.  That  is, 
cosh  y(*,  ♦ * - b) 


B(h,z)  ^ 


J- 

cosh  Yfj 


2 cosh  Yt  sinh  y(i^j  ■ b) 
sinh  Y<lj 


(58) 


Furthermore, 


0 8Blh..l  • * 

Y 3z  sinh  y^i 


h) 


2 


sinh  Y*  sinh  y(<^j 
sinh  Yflj  tanK  yAj 


(59) 


VERTICAL  MAGNETIC  DIPOLE  (VMD)  ANTENNA 

By  employing  Maxwell's  equations  and  the  reciprocity  theorem,  we  see 
that  the  quasi-near  range  V'MD  antenna  field-component  expressions,  which  are 
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valid  within  the  conducting  slab,  may  easily  be  determined  from  the  HED 
antenna  expressions.  They  are 


Q A(h,z)  . 


H Q B(h.z)  , 

O “ 


..  9IA  ^2  , ,, 

”2  TT  Q A(h,2)  , 

2irY  p 


where  A(h,2)  » A(2,h),  with  z and  h interchanged.  That  is. 


A(h,z)  A. 


2 cosh  Yz  cosh(tj  - h)  sinh  yC^^  z - h) 


cosh  Y^i 


sinh  yl^  coth  y^ j 


niSCUSSION 


It  would  he  of  interest  to  compare  the  results  derived  in  this  report 
with  some  known  results.  Figures  2 through  6 show  comparisons  of  modified 
image  theory,  quasi-near  (asymptotic)  theory,  and  Weaver's  numerical  integra- 
tion results"  for  the  electric  nnd  magnetic  fields  at  the  surface  (z  » 0) 
of  a one-skin-depth-thick  (Ij  «=  6)  conducting  slab  produced  by  an  HED  located 
in  the  middle  of  the  slab  (h/6  « 0.5).  The  normalized  amplitude  of  each 
component  (E'  or  H')  is  given  by 


It  * 

I cos  ♦, 


and  E' 


— ~ ts 

/sin  ^ \ * 

Wos 


Two  values  of  a and  b are  considered  for  the  modified  image  theory  plots. 
The  a « 0.4  and  b « 0.96  results  should  be  valid  in  situations  close  to  the 
source,  and  the  a ■ 0.96  and  b » 0.4  results  should  be  valid  at  further 
distances.  In  each  of  these  cases,  only  five  terms  of  the  infinite  summation 
were  needed  for  1 percent  accuracy. 

As  Weaver  has  indicated,®  this  particular  model  possesses  s)Tnmetry 
about  the  plane  z ■ h.  When  z - h,  all  components  that  vary  as  B(z,h)  equal 
zero.  Furthermore,  if  z ■ 0 or  ti,  all  components  that  vary  as  A(z,h)  are 
equiil,  whereas  all  components  that  vary  as  B(z,h)  are  equal  and  opposite. 
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Figures  2 and  3 show  horizontal  electric  field  coaparisons  of  the  modi- 
fied image  theory,  asyaptotic  theory,  and  numerical  integration  results.  We 
can  see  from  figure  2 that,  for  the  Ed  component,  the  modified  image  theory 
a B 0.4  and  b ■ 0.96  curve  agrees  well  with  the  numerical  integration  results 
throughout  the  range  of  p/d  considered  (0  ^ p/d  < 3).  However,  we  observe  { 

that  beyond  approximately  1.25  skin  depths7  the  Tsinple  form)  asymptotic 
theory  agrees  more  closely. 

A ^ 

We  can  see  from  figure  3 that,  for  the  E4  component,  the  modified  image 
theory  a ■ 0.4  and  b ■ 0.96  curve  agrees  well  with  the  numerical  integration 
results  for  0 ^ p/d  ^ 0.75,  while  the  a » 0.96  and  b « 0.4  curve  more  closely 
agrees  for  0.75  < p/d  ^ 3.  Beyond  approximately  2 skin  depths,  the  asymptotic 
theory  provides  the  better  fit  to  the  numerical  integration  data. 

Figures  4,  5,  and  6 show  the  magnetic  field  comparisons  of  the  modified 
image  theory,  asynptotic  theop^,  and  numerical  integration  results.  We  see 
from  figure  4 that,  for  the  Up  component,  the  a ■ 0.4  and  b ■ 0.96  curve 
agrees  very  closely  with  the  numerical  integration  results  for  p/6  ^ 0.5,  but 
only  moderately  so  for  p/d  > 0.5.  The  a ■ 0.96  and  b » 0.4  curve  is  in  good 
agreement  beyond  1.5  skin  depths,  while  the  asymptotic  theory  curve  provides 
the  best  fit  beyond  2.5  skin  depths. 

In  figure  5,  we  can  see  that,  for  the  H4  component,  the  modified  image 
theory  a ■ 0.4  and  b » 0.96  curve  agrees  very  closely  with  the  numerical 
integration  results  to  p/d  'v-  2.  Beyond  that,  the  asymptotic  theory  curve 
provides  the  best  fit. 

In  figure  6,  we  can  see  that,  for  the  Hz  component,  the  modified  image 
theory,  a « 0.4  and  b ■ 0.96  curve  agrees  very  well  with  the  numerical  inte- 
gration results  for  p/d  < 1.  Beyond  p/d  ■ 1,  the  a ■ 0.96  and  b ■ 0.4  curve 
agrees  more  closely.  For  p/d  > 2,  the  asymptotic  theor>'  curve  also  agrees 
well  with  the  numerical  integration  results. 

Thus,  it  appears  that  the  (simple  form)  asymptotic  theory  will  provide 
results  of  sufficient  accuracy  when  the  measurement  distance  is  greater  than 
2 skin  depths  and  greater  than  twice  the  slab  depth  (i.e.,  p/d  > 2 and 

p/tj  > 2). 

CONCLUSIONS 


Approximate  expressions  for  both  the  general  ac  quasi-static  and  quasi- 
near  fields  produced  by  electric  and  magnetic  dipole  antennas  located  within 
a conducting  slab  have  been  derived  by  employing  finitely  conducting  earth- 
image  theory  techniques  and  by  applying  the  quasi-near  approximation  to  the 
basic  Sommerfeld  integrals. 

We  have  demonstrated  that  the  resultant  approximations  very  closely 
agree  with  previously  derived  numerical  integration  results.  In  particular, 
it  appears  that  the  (simple  form)  asymptotic  theory  will  provide  sufficiently 
accurate  results  when  the  measurement  distance  is  greater  than  2 skin  depths 
and  greater  than  twice  the  slab  depth. 
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Although  displcceMiit  currents  In  the  conducting  slab  have  been  ignored 
in  the  analysis,  they  can  be  included  by  siB|>ly  replacing  o with  o ♦ iwe  in 
the  field  strength  equations,  providing  |y*|  >>  |yq*|.  The  resultant  expres- 
sions are  applicable  to  short  range  electromagnetic  propagation  in  a shallow 
sea. 
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